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We report the cooling of electrons in nanoelectronic Coulomb blockade thermometers 
below 4mK. Above 7 mK the devices are in good thermal contact with the environ¬ 
ment, well isolated from electrical noise, and not susceptible to self-heating. This is 
attributed to an optimised design that incorporates cooling hns with a high electron- 
phonon coupling and on-chip electronic hlters, combined with a low-noise electronic 
measurement setup. Below 7 mK the electron temperature is seen to diverge from 
the ambient temperature. By immersing a Coulomb Blockade Thermometer in the 
^He/^He refrigerant of a dilution refrigerator, we measure a lowest electron temper¬ 
ature of 3.7 mK. 


Electronic mail: 
Electronic mail: 
Electronic mail: 


j. prance@lancaster .ac.uk 
leif. roschier@iki. fi 
mika.prunnila@vtt.fi 


1 



Cooling nanoelectronic structures to millikelvin temperatures presents extreme challenges 
in maintaining thermal contact between the electrons in the device and an external cold 
bath^. It is typically found that the electrons in a nanoelectronic device cooled even to 
~ 10 mK are signihcantly overheated. Understanding how to obtain and measure elec¬ 
tron temperatures approaching 1 mK will open a new regime for studying nanoelectronics, 
and pave the way towards pioneering sub-millikelvin techniques^. This would beneht, for 
example, investigations of the fractional quantum Hall effect in 2D electron gases^’"^, and 
solid-state quantum technologies including superconducting and semiconducting qubits. To 
access these temperatures, one must minimise parasitic heating and internal Joule heat¬ 
ing, maximise the coupling to cold contact wires and phonons in the host lattice, all the 
while overcoming the decrease in electron-phonon coupling and electrical heat conduction 
as temperatures drop. 

Here we study Coulomb Blockade Thermometers (CBTs) that have been designed to 
operate signihcantly below 10 mK and demonstrate cooling of electrons in a nanoelectronic 
device below 4mK. This was achieved using CBTs as a diagnostic tool to quantify and 
optimise the thermal environment of the electrons in the device. A CBT consists of an array 
of Coulomb-blockaded metallic islands connected by tunnel junctions whose conductance 
is temperature dependent^’®. CBTs typically function over a decade of temperature and 
have previously been demonstrated to work at temperatures as low as 7.5 mK^. Perhaps 
most importantly, they can be viewed as a primary thermometer of their internal electron 
temperature. 

The structure of the CBTs studied here is shown in Fig. 1. Devices are fabricated using an 
ex situ tunnel junction process® that provides excellent tunnel junction uniformity®, and has 
also been used to fabricate superconducting qubits^®. In this work the CBTs are made from 
A1 hlms with a thickness 250 nm and an isolating layer of 250 nm SiOx deposited by PECVD. 
The tunnel junctions have a nominal diameter 0.8 pm and a resistivity lOkDpm^. The 
substrate is undoped Si with 300 nm thermal oxide on the surface. 

Efficient thermalisation of electrons and phonons in the metallic islands of the CBT is 
critical for reaching low electron temperatures^^’^^. The electron-phonon heat flow Pep is 
described by the material-dependent electron-phonon coupling constant S and the volume 
of the metallic island D, 

ftp = ESI (p - T=) (1) 
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FIG. 1. Optical micrograph of the CBT with equivalent circuit diagram (left) and schematic 
cross-section of the structure (right). The CBT is formed of 32 x 20 metallic islands of capacitance 
Cy, connected in an array by tunnel junctions of resistance Rt, as shown in the circuit diagram. 
Connection to the array is made via on-chip RC filters comprising a meandering electrode sand¬ 
wiched between large-area grounded metal films, separated by 250 nm SiOx. Each filter has a 
distributed resistance R ~ 500 0 and capacitance C ~ 10 pF. The schematic cross-section (right) 
shows one tunnel junction connecting two islands, with Au thermalisation blocks on top. 

where Tg is the electron temperature and Tp is the phonon temperature^^. To minimise 
Tg, the island volume should be large and the material chosen to maximise S. The island 
thickness achievable with the ex situ tunnel junction process, or other deposition techniques 
used for tunnel junction devices, is typically up to 1 pm. Thicker films suffer from stress 
build-up, causing poor adhesion between the film and the substrate. This is a severe prob¬ 
lem at mK temperatures where poor adhesion can lead to poor thermalisation and even 
mechanical failure due to thermal motion during cool-down. We avoid these problems by 
using a combination of the ex situ process followed by masked electroplating of Au on top 
of the CBT islands^'^, which we refer to as thermalisation blocks. Electroplating can pro¬ 
duce ~ 10 pm thick, low stress Elms, and here we choose a nominal thickness of 5 pm for the 
thermalisation blocks. The effective electron-phonon coupling in these islands, with nominal 
volume 5 x 205 x 38.5 pm^ and a high coupling constant^^ in Au S = 2.4 x 10® W m“^, is 
estimated to be more than two orders of magnitude larger than in previous CBTs fabricated 
using the ex situ junction process^®. 

In addition to the thermalisation of the CBT itself it is important to thermalise the 
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incoming leads through robust thermal anchoring and heavy electromagnetic hltering^^. We 
improve the chain of thermalisation and hltering by including on-chip resistive meander 
structures in line with all electrical contacts. These form a distributed resistive-capacitive 
chain with a cut-off frequency of ~ 40 MHz. Similar hlters based on a large area capacitor 
and tunnel junctions in series have previously been incorporated in a CBT^®. 

Figure 2 shows the behaviour of a CBT, fabricated using the process described above, 
focusing on temperatures between 7 mK and 80 mK. The sensor was measured in a com¬ 
mercial cryogen-free dilution refrigerator^® with a base temperature ~ 7 mK and later in a 
custom dilution refrigerator manufactured at Lancaster University^® with a base tempera¬ 
ture ~ 2.5 mK (see below). The conductance of the CBT was measured in a current-driven 
four-wire conhguration, with the drive current and voltage amplihcation provided by an 
Aivon PAIO ampliher. A small AC excitation (typically 5pA < Jag < 100 pA) was added 
to the DC bias Ibc, allowing the differential conductance G to be measured with a lock-in 
ampliher. 

As shown in Fig. 2(a), the CBT conductance dips around zero bias, and the dip becomes 
deeper and narrower at lower temperatures. Its full-width at half maximum is related 
to temperature by V 1/2 ~ bASdNksT/e, where N is the number of tunnel junctions in 
series^. This does not account for self-heating in the sensor and so is only applicable when 
T = Te = Tp. A more practical parameter to determine temperature is the zero-bias 
conductance Go, which has an approximate analytic relation to temperature^® 

Go ^ Gt (1 - 1/Qun - l/60u% + 1/630m^) (2) 

where un = Ec/^bT is dimensionless inverse temperature, Ec = [(iV — l)/iV]e^/Gs is the 
charging energy of the system, Gs is the total capacitance of each island, and Gt is the 
asymptotic conductance. When Un < 2.5 the temperature measurement error is < 2.5%^^. 
Thus, if Gs and Gt are known it is possible to determine T by measuring only Go. The 
most complete method to determine temperature is using a full tunnelling model to calculate 
G(Vbc) numerically^. We use this last approach to hnd Gs and Gt for the device in order 
to then determine the temperature. 

Numerical calculations of conductance are made using an algorithm derived from the 
free, open-source library pyCBT^^. In addition, we account for overheating in the sensor by 
predicting the electron temperature Tg in the islands using a model for the heat flow into 
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FIG. 2. CBT behaviour between 80 mK and 7 mK combining measurements made in two dilution 
refrigerators, (a) CBT conductance G versus measured bias voltage Vdc at four temperatures. 
Symbols show measured values and lines show best fits to the calculated ideal conductance (see 
text for details of fitting). The warmest measurements (crosses, circles and squares) are fitted 
simultaneously to calibrate the sensor, giving Cs = 236.6 fF and Rt = 22.42 kfl. The coldest 
measurement (triangles) was htted using this calibration. The minimum electron temperatures 
are in close agreement with the refrigerator temperature measured by the Ru 02 thermome¬ 
ter: 59.9 mK, 40.1 mK, 20.0 mK, and 7.2 mK respectively, (b) CBT electron temperature Tg vs 
refrigerator temperature Tmxc- Symbols show from fits to conductance dips measured in the 
cryogen-free refrigerator (circles) and the custom refrigerator (triangles). The solid line shows Tg 
determined by monitoring the conductance at I^DC = 0 as the cryo-free fridge cooled over 35 min, 
showing that the CBT has a stronger thermal link to the refrigerator than the Ru02 thermometer, 
leading to the thermal lag (Tmxc > Re) during this time. 
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each island, 


P = ^ - Sfi (T^ - T^) + Po (3) 

ilj’ 

where the hrst term is Joule heating at tunnel junctions of resistance P^, the second term is 
heat flow to phonons, and Pq accounts for parasitic heating. For a given phonon temperature, 
the minimum electron temperature is found at Voc = 0. If Pq is small, then ^ Tp. 

Figure 2(a) shows the result of htting the calculated G(Vdc) simultaneously to three 
measurements made between 20 mK and 60 mK. The £t parameters are Gt, C*s, and Tp for 
each measurement. The fit was found to be insensitive to the value of Pq, and hence parasitic 
heating is assumed to be negligible at these temperatures. Having calibrated the CBT, the 
fitted P's was used to fit further measurements with Tp and Gt as the free parameters. An 
example is given in Fig. 2(a), where the electron temperature is found to be 7.3mK at a 
refrigerator temperature of Tmxc = 7.2 mK, here measured by a calibrated Ru02 resistor^^. 

Figure 2(b) shows that the electron temperature measured by the CBT is in close agree¬ 
ment with the refrigerator temperature T^xc in both refrigerators between 7 mK and 80 mK. 
In the custom refrigerator, Tmxc is measured using a conventional Vibrating Wire Resonator 
viscometer (VWR) immersed in the saturated dilute phase of the ^He/^He refrigerant in 
the mixing chamber^^’^^. For this set of measurements, in both refrigerators the CBT is in 
vacuum and housed in a gold-plated copper package (Aivon Oy SH-1) that is attached to the 
mixing chamber plate. The package includes RC filters with a cut-off frequency ~ 300 kHz. 
Electrical contacts to the CBT are thermalised in additional cold RG filters potted in Ec- 
cosorb CR 124 (Aivon Oy ‘Therma’), which are also attached to the mixing chamber plate. 
The filters have a cut-off frequency ~ 15 kHz. 

Subsequently we investigated the behaviour of the CBT below 7 mK in the custom dilution 
refrigerator. Eigure 3(a) shows the electron temperature of the sensor, when in vacuum, 
gradually cooling below 4 mK while the refrigerator was held at Tmxc ~ 2.8 mK. Here Tg was 
determined from the value of Go, which was found by measuring G over a small range of IJdc 
(~ 30 pV) close to Vdc = 0. We observe an extremely long equilibration time (over 3 days) 
but a rapid cooling of the CBT following a heating event (a rehll of the liquid helium bath 
that briefly heated the CBT above 5.5 mK and the fridge above 3.5 mK). This suggests that 
the thermal contact between the CBT and the refrigerator is relatively strong, and that its 
slow cooling is due to heat leak from an external warm object. 

A second CBT was immersed in the ^He/^He refrigerant of the custom dilution refrigerator 
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FIG. 3. Thermalisation of two CBTs as a function of time at a refrigerator temperature of 
Tuoic < 2.8 mK. (a) shows cooling of one CBT in vacuum (circles) and one immersed in the 
^He/^He refrigerant of the dilution refrigerator (triangles). In both cases, the CBTs are cooling 
after being warmed above 10 mK by temporarily increasing the refrigerant temperature. The CBT 
in vacuum is extremely slow to thermalise. By comparison, the CBT immersed in ^He/^He ther- 
malises significantly faster. The inset shows cooling of the immersed CBT after it has been heated 
by a large DC drive current (50nA, 40nA and 30nA for run 1, 2 and 3 respectively). Fitting to 
an exponential decay (solid line) yields a time constant of 570 s and a saturation temperature of 
3.8mK.(b) shows schematic of the immersion cell. 

to improve thermal coupling and better isolate external sources of heating. A schematic 
of the immersion cell is shown in the Fig. 3(b). Sintered silver blocks increase the thermal 
contact between the refrigerant and the CBT and contact wires. Several sinters was attached 
to the sensor package, to each of the four measurement wires, and a grounding wire for the 
package and the RC filters. The immersed CBT was found to equilibrate much faster, as 
shown in Fig. 3(a), reaching Te = 3.8 mK at T^xc = 2.7mK. 

In order to study the CBT at |Vdc| 0, the time needed for the sensor to reach thermal 
equilibrium after a change of Joule heating needs to be known. The inset to Fig. 3(a) shows 
the relaxation of Tg after the CBT has been heated by a large drive current for long enough to 
reach thermal equilibrium (> 30 min). The subsequent value of Tg is measured by scanning 
close to Voc = 0, where Joule heating should be negligible. The relaxation of Tg is found to 
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FIG. 4. Characteristics of a CBT immersed in ^He/^He refrigerant, (a) Fitting to the warmest three 
measurements gives Gs = 209.5 fF and Rt = 23.21 kfl. The fitted minimum electron temperatures 
for the warmest three curves are in reasonable agreement with the refrigerator temperature as 
measured by the VWR thermometer: 29.4 mK, 19.0 mK, 10.5 mK respectively, (b) shows measured 
electron temperature in the CBT as the refrigerator cooled steadily from 10 mK to 2.7 mK over a 
period of 12 hours. The solid line shows a fit of the form T* = T^xc + c, yielding an exponent 
X = 2.7. The dot-dashed line shows a best fit of = TJ^xc + c. The inset to (b) shows the 
measured Te as the refrigerator was temporarily cooled in single-shot mode to 2.2 mK, reaching a 
lowest Te below 3.7 mK. 

have a time constant of 570 s. 

Fignre 4(a) shows the calibration of the immersed sensor. The three warmest measnre- 
ments were htted simnltaneously to determine Cs = 209.5 fF and Rt = 23.21 kC. The htted 
temperatures agree with T^xc to within 6%. Given the agreement between the htted Tg™™ 
and Tmxc, we can assume that parasitic heating is still negligible down to 10 mK. 





























The coldest measurement in Fig. 4(a) was fitted using the above values, yielding a mini¬ 
mum electron temperature of 3.9 mK. This measurement was made over a period of 7 hours 
to ensure that the CBT was in thermal equilibrium at each value of Vdc- At these tem¬ 
peratures the parasitic heating of the CBT is now signihcant and does not match the 
refrigerator temperature of T^xc = 2.7 mK. In order to £t this conductance dip the ther¬ 
mal model, Eq. 3, was used with Tp = T^^c, and with the parasitic heating Pq and the 
electron-phonon coupling Shi as free parameters. 

The overheating of the sensor at Vdc = 0 constrains the product however, the 

parasitic heating is not large enough to reliably separate Pq and SV2 in the fit. Qualitatively, 
the fits suggest that Pq > 300 aW per island and EV2 is at least four times larger than 
expected from the nominal size of the thermalisation blocks and the literature value of S for 
Au^^. It was not possible to determine an upper bound on Pq. It is worth noting that the 
power required to measure the CBT conductance (< 1 aW per island due to Joule heating 
from Jac) is much lower than our estimate of Pq. As such, we believe that CBTs of this 
type can be operated at still lower temperatures by reducing parasitic heating. 

Figure 4(b) shows how the CBT electron temperature diverges from the refrigerator 
temperature below ~ 7 mK. Here the value of Te is found by measuring Gq close to Vdc = 0 
and so Joule heating can be neglected. The functional form of Te vs T^xc should have 
the same temperature dependence as the dominant thermalisation mechanism, i.e. T® for 
electron-phonon coupling. However, other power laws have been observed^. Here we find 
that the best fit of + c gives x = 2.7 and a saturated Te of = 3.4 mK. The 

fitted exponent x cannot be confirmed by measurements of the conductance dip because the 
overheating is still relatively weak even at the lowest temperatures. We find that a thermal 
model with a T^ thermalisation term fits the data equally well as a model using T®. In order 
to understand the thermalisation mechanism in more detail, this sensor would need to be 
cooled closer to 1 mK. 

In conclusion, CBTs of the structure described here have been shown to operate as reliable 
primary thermometers of electron temperature down to 3.7mK. The large thermalisation 
blocks incorporated in the device, and a relatively low level of parasitic heating, ensure 
that the electron subsystem in the sensor is well coupled to the phonon subsystem down to 
^ 7 mK. An immersion cell was shown to improve thermal coupling between a CBT and a 
dilution refrigerator. This allowed the onset of overheating to be observed below 7 mK, and 
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while the presence of overheating could be seen, the effect was sufficiently weak that the 
sensor will need to be cooled further to fully characterise the thermalisation mechanisms. 
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